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ABSTRACT 

The continuity and heat conduction equations for 
electrons, ions and neutral species are solved simultaneously 
from the viewpoint of studying the F-region magnetic storm 
behavior. It is shown that many of the observed changes 
in the F-region during the main phase of a magnetic storm 
can be successfully explained by assuming a decrease in oxygen 
concentration at the turbopause level. The specific changes 
in neutral composition, electron density, and electron , ion , 
and neutral temperatures during magnetic storms are discussed 
for different solar conditions. 


*To be published in PLANETARY AND SPACE SCIENCE 



Introduction 


The storm time morphology of the bottomside F-region has 

been extensively studied during the past three decades [see 

Obayashi (1964) for detailed references]. The detailed behavior 

of the storm time F-region cannot be described in a simple way, 

since it depends upon the phase and time of the storm as well 

as the location of the observing station. However, the nearly 

universal features of a magnetic storm are a worldwide decrease 

in foF2 (up to 50 percent) and the increase in the height of 

the electron density maximum (up to 100 km) . The electron 

density distribution in the topside region during disturbed 

conditions has been studied only recently [Evans (1965), Reddy 

et.al (1967), Hibberd and Ross (1967), Rao (1968), and Bauer 

and Krishnamurthy (1968)]. Here again, though the detailed 

behavior depends on the height, location and phase of the storm, 

in general the electron density rises at all heights from about 

100 km above the F2 peak. In addition, the total electron 

content decreases in a manner similar to n F2. 

m 

Since the satellite drag data has revealed a worldwide 
enhancement in the neutral density and neutral temperature 
(the latter derived from the density data by comparison with 
an atmospheric model) during magnetic disturbances (Jacchia 
et.al. 1967) it was considered likely that the depletion in 
a m F2 may be due to enhancement in the upper atmospheric temperature 
The actual heating could be caused by one of the several mechanisms 
hydromagnetic waves (Dessler 1965) , heat conduction waves 
(Volland 1967) and Joule dissipation of the ring current (Cummings 
and Dessler , 1967) . 

Garriot and Rishbeth (1963) and Thomas (1966) argued on 
the basis of solutions to the electron density continuity 
equation that the decrease in n F2 caused by atmospheric heating 
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cannot be more than 25 percent and hence not sufficient to 
account for the large decrease in density actually observed,, 
Several theories exist for the F-region disturbances. The 
electrodynamic drift-theory, (Martyn , 1954) , the increase in 
dissociative recombination rate (Seaton 1956) , the flow of 
ionization between ionosphere and magnetosphere through 
hydromagnet ic pumping (Piddington 1964) and the loss of thermal 
plasma through inflation of the magnetosphere and through 
local "energization" (Bauer and Krishnamurthy , 1968), Many 
of these ideas have certain merits since they succeed in 
explaining the general storm time behavior of the F-region 
to varying degrees of approximation. The difficulty comes, 
however, when an attempt is made to account for the temperature 
and density variations of both the ionized and neutral consti- 
tuent throughout the altitude range of the F-region, Moreover, 
many of these ideas are qualitative and it is difficult to 
relate them to the observational data in a quantitative manner. 

We feel, however, that the temperature effects during the 
disturbed condition are quite important in spite of the dif- 
ficulties encountered by Garriot and Rishbeth (1965) and Thomas 
(1966), Ira this paper, we shall attempt to examine this aspect 
of the F-region magnetic storm. 

The Effect of Changes in the Lower Boundary 

In the accompanying paper [Herman and Chandra 1968; to be 
referred to as paper l] the effects of' solar EUV on density and 
temperature of both the neutral and ionized constituents were 
discussed. It was shown that in solving the equations of 
continuity and heat transport, the boundary values of the 
neutral constituent cannot be assumed invariant throughout the 
solar cycle. In discussing the magnetic storm problem, a similar 
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situation arises. It is reasonable to assume that the solar 
EUV flux maintains its pre-storm values, but there is no reason 
to believe that the lower boundary densities remain invariant 
during the disturbed period. The neutral composition , part icularly 
atomic oxygen, in the mesosphere and the lower thermosphere are 
strongly controlled by eddy diffusion, transport, and chemical 
processes [Colegrove et . al. 1966, Shimazaki 1967, Maeda and 
Aikin 1968, Hesstvedt 1968], In the high latitudes these 
processes can be modified by the precipitation of energetic 
electrons and the heat dissipation due to ring current. Hines 
(1965) believes that some of the energy deposited in the auroral 
region during magnetic storms can be transported to the lower 
latitudes by means of gravity waves. While it is difficult to 
estimate the composition changes resulting from these mechanisms, 
it is reasonable to postulate that the composition in the lower 
atmosphere does change during magnetic storms particularly in 
the middle and high latitudes. Evidence of such changes in 
composition can be inferred from the enhancement in the intensity 
of the airglow during magnetic storms [Truttse, 1968; Weill 
and Glaume, 1967]. It will be shown in the following that many 
of the important features of the F-region magnetic storm can 
be described from the simultaneous solutions of the equations 
of continuity and heat transport for both the neutral and ionized 
constituents by making suitable changes in the lower boundary. 

The Equations of Continuity and Heat Conduction 

The equations of continuity and heat transport and their 
method of solution have been discussed in paper I. In the 
following, only an outline is presented for ready reference. 

Equation of continuity for electron density: 
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Heat conduction equation for electrons: 
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Heat conduction equation for ions: 
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Heat conduction for neutral constituents: 
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The various neutral constituents (N 2 , Og, 0, H e , and H) are 
assumed to be in diffusive equilibrium 
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Discussion 

The definition of the terms in equations 1-5 and their 
numerical values are exactly the same as in p^per I as are their 
method of solution. Here we have characterized the storm time 
conditions by a decrease in the oxygen concentration at the turbo- 
pause level as compared to the pret-storm value. As in paper I, 
this level has been assumed to be 100 km. It may be pointed out 
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that in many of the atmospheric models, the turbopause level 
is generally assumed to be between 120-130 km. Our present 
choice has been guided by the recent experimental results 
(Shaeff er 1967; Hartman et.at. 1967) which show this level to 
be less than 110 km. The general conclusion of this paper is 
not dependent on the exact knowledge of the reference boundary 
values. 

The numerical results of the calculations are shown in 
Figure 1-3 which show the steady state height profiles of the 
densities and temperatures for both the neutral and ionized 
constituents corresponding to the various phases of the solar 
cycle. We have chosen three values of EUV fluxes, 1.00 HF, 

1.30 HF and 1.60 HF, as representative of different phases of 
solar activity (HF denotes the flux values in the spectral 
range of 10-1027A as given by Hinter-egger et.al. (1965)). 

All the calculations relate to the magnetic dip angle 70° 
and solar zenith angle 39°, the latter corresponding to approxi- 
mately equinoctial noon condition at mid-latitudes. 

Figure (1-a) shows the density profiles for [Ng] ? [0g] 
and [ 0 ] corresponding to 1.30 HF. The boundary values assumed 
for the five neutral constituents for the quiet and disturbed 
condition are as given in the following table. 


N r 


0 


2 


H 


H 


TABLE 1 
QUIET 

4.65 x 10^ 2 /cm 3 
1.25 x 10 12 /cm 3 
1.00 x 10 12 /cm 3 
2.50 x 10 8 /cm 3 
1.00 x 10 3 /cm 3 


DISTURBED 
4.65 x 10 12 /cm 3 
1.25 x 10 12 /cm 3 

5.00 x 10 1 ' L /cm 3 
2.50 x 10 8 /cm 3 

1.00 x 10 3 /cm 3 
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It may be noted that the only difference assumed in the quiet 
and disturbed condition is the reduction in oxygen concentration 
by 50 percent compared to its pre-storm value at the reference 
level* The values assumed for the two cases are quite arbitrary 
and have been adopted for illustration only. The density and 
temperature profiles as a result of this decrease are shown in 
Figures 1-a, l*=b„ The main points of these Figures can be 
summarized in the following: 

1» The isothermal neutral temperature increases from about 

900°K to 1100°K - an increase of about 200°K (Figure 1-b) . 

2s There is a general increase in [Ng] and [Og] values 
throughout the altitude range corresponding to the 
increased neutral temperature (Figure 1-2). The helium 
and hydrogen densities show a similar increase (not 
shown here for the sake of clarity) . 

3« The oxygen concentration is less than its pre-storm value 
up to about 300 km though the difference between the 
two values [°]q - [0]^ continuously decreases with 
altitudes At about 300 km the two profiles cross over 
and from there on [0]^ continues to be greater than [0 ]q» 

4. Both the electron and ion temperatures show an enhancement 
by a few hundred degrees compared to the quiet condition. 
The peak electron temperatures T e maxQ and T e max D are 
2170 a»d 2479 ° k at 180 and 200 km respectively. 

5 0 The electron density n eD is less than n 0 Q up to about 

350 km. The two distributions cross over at this altitude. 

Above 350 km n eD continues to be higher than n 0 Q. The 

maximum decrease in n e p is near the height of the maximum. 

The maximum electron densities n and n. n are 

emaxQ emaxD 
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respectively 13,5 x 10^ and 6.8 x 10^ elect. /citf* at 

This reflects a decrease of about 


260 km and 310 km. 

50 percent in n 

The total electron content n 


„ and increase in h „ by 50 km, 

emax max 


and n. n are 

19 -9 tQi 9 tD 

23.2 x 10 and 19.4 x 10 x /cm“2. This shows a 

decrease of about 16 percent. 


Figures 2a and 2b show the results corresponding to 1.00 HF 

which corresponds to low solar activity. The important features 

of these figures are similar to the ones discussed in the 

preceeding section though the changes in the disturbed and quiet 

conditions are not quite the same. The important differences 

are in T g and n g profiles. The electron temperature shows a 

transition from a peaked profile; to a monotonically increasing 

profile (Figure 2b). The increase in T g is almost by a factor 

of 2 from about 200 km onwards. The effect of this considerable 

increase in T is reflected in the n profile which shows a rather 
e e 

large gradient above the F^ peak. This abnormal increase in 
T e should not be interpreted as typical of a mid»latitude storm 
during low solar activity. The increase is usually associated 
with the change in the characteristic of the T q profile as was 
discussed in paper I. In fact, on either side of the transition 
region the percentage variation in T 0 is rather small. 

Figures 3a, and 3b show the density and temperature profiles 
corresponding to HF = 1.6. The main features are essentially 
the same as those shown in Figures la and lb. The following 
table summarizes the overall results of Figure 1, 2 and 3. 
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Table 2 shows some numerical examples of the effects of 
decreasing oxygen concentration on the various ionospheric 
parameters. These changes are consistent with the commonly 
known changes in the F-region during magnetic storms in the 
middle and high magnetic latitudes. Thus, our calculations 
show reasonable agreement with the observed enhancement of 
electron density in the topside region of the order of 100 
percent [King et.al., (1967) Baiier 1 and Krishnamurthy (1968) ] and the 
depression in n F„ up to 50 percent (Obayashi 1964) . The 
increase in the height of the F 2 maximum, electron, ion, and 
neutral temperatures and neutral density, and the decrease in 
total electron content as shown in Table 2 [and Figure 1-3], 
are in general agreement with the observational data [Evans 
1965, Jacchia et.al., Hibberd and Ross, 1967], We cannot, of 
course, explain the complex changes in the F-region through 
the various phases of magnetic storms based on steady state 
solutions presented in this paper. For a complete understanding 
of the problem it will be necessary to consider the solutions 
of the time dependent equations along with the time dependent 
boundary values appropriate for storm time conditions. 

In our calculation we only decreased the oxygen concentration 
and the other changes followed accordingly. The result would 
have been the same if instead of decreasing [0] at the turbo- 
pause level, we would have lowered the level itself. Further, 
we did not change the relative values of [0] with the change in 
solar activity, for the sake of simplicity. The main purpose of 
this paper is to illustrate the result of this change rather than 
to obtain a quantitative agreement between the theoretical and 
observational results. 
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Conclusion 

In this paper we started with the assumption that the 
turbopause level of oxygen is decreased during the main phase 
of the magnetic storm in the middle and high magnetic latitudes 
resulting in the decrease of [0] in the lower thermosphere. 

The general behavior of the F-region during magnetic storms can 
be satisfactorily explained with this assumption. Consequently 
the decrease in foF2 and the increase in the neutral temperature 
and density are the results of redistribution of [0]. 

The physical implications of the changes in oxygen concentra- 
tion is easy to understand if it is realized that the atmospheric 
neutral temperature is essentially controlled by the relative 
magnitude of heat production and heat loss in the lower thermosphere. 
The thermal conduction helps only to establish isothermal tempera- 
tur<es at high altitudes. The electron density distribution is 
controlled by the photo-chemical processes in the lower F-region 
and photo-diffusive processes in the upper F-region. Since 
oxygen provides an important heat sink for the neutral constituents 
in the lower thermosphere, a reduction in its value results in 
the reduction of the loss term in the neutral temperature equation. 
The contribution to the heat production term is also reduced, 
though not to the same extent as the loss term. Consequently, 
there is a net increase in the total heat input to the neutral 
gas, and hence, an increase in neutral temperature. The various 
neutral constituents, which are assumed to be in diffusive 
equilibrium, are also increased corresponding to the new higher 
temperature. In the case of [0], this increase is naturally 
reflected only in the higher altitudes [Figures 1-a, 2-a, and 3-a] . 
The effect of this redistribution is to increase the electron 
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loss term [due to its direct dependence on [Ng] and [0^]] and 
to decrease the electron production term [due its direct 
dependence on [0]] in the bottomside F~region. The net effect 
then is an overall reduction in electron density. In the 
topside region the situation is somewhat different. Here the 
electron production term is increased and the diffusion coefficient 
is decreased. Both these changes have a net effect of increasing 
the electron density. 

Unfortunately, the observational data on the neutral 
composition in the lower atmosphere during magnetic storms are 
almost non-existent. Recently, Reber (private communication) 
from his mass spectrometer experiment aboard Explorer XXXII, 
found some evidence of composition changes during magnetic 
storms. His observations revealed a decrease in [0] by 45 
percent and increase in [Ng] by a factor of 2 at 285 km and 
in the latitude (geographic) range of 36°-46°. This observa- 
tion is in general agreement with our result. Clearly many more 
observations of this kind are needed particularly in the lower 
thermosphere and mesosphere for a direct experimental verification 
of our hypothesis. It is evident, however, from the present 
study that the changes in the F~region are closely related to 
the composition changes in the lower atmosphere. An understanding 
of the photo-chemical and dynamical processes effecting this 
region both during quiet and disturbed conditions is therefore 
of considerable importance in understanding the storm time behavior 
of the F-region. 
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The altitude variations of Ng, Qg and 0 for 
quiet and disturbed conditions corresponding 
to HF - 1.3. The oxygen density variations 
between 100 to 300 km have been redrawn in 
the upper right hand corner for clarity , 

The altitude profiles of n , T , T. and T 

0 6 x n 

for quiet and disturbed conditions corresponding 
to HF = 1.30. 

The altitude variations of Ng, 0g and 0 for 
quiet and disturbed conditions corresponding 
to HF = 1.0. The oxygen density variations 
between 100 to 300 km have been redrawn in 
the upper right hand corner for clarity . 

The altitude profiles of n , T , T. and T 
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for quiet and disturbed conditions corresponding 
to HF = 1.0. 

The altitude variations of Ng, 0^ and 0 for 
quiet and disturbed conditions corresponding 
to HF - 1.6. The oxygen density variations 
between 100 to 300 km have been redrawn in 
the upper right hand corner for clarity. 

The altitude profiles of n , T , T. and T 
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for quiet and disturbed conditions corresponding 
to HF — 1.6. 
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Figure 2-b . 
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